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Abstract: We develop a set of controlled, analytic approximations to study the effects 
of bubble collisions on cosmology. We expand the initial perturbation to the inflaton field 
caused by the collision in a general power series, and determine its time evolution during 
inflation in terms of the coefficients in the expansion. In models where the observer's 
bubble undergoes sufficient slow-roll inflation to solve the flatness problem, in the thin 
wall limit only one coefficient in the expansion is relevant to observational cosmology, 
allowing nearly model-independent predictions. We discuss two approaches to determining 
the initial perturbation to the inflaton and the implications for the sign of the effect (a hot 
or cold spot on the Cosmic Microwave Background temperature map). Lastly, we analyze 
the effects of collisions with thick- wall bubbles, i.e. away from the thin- wall limit. 
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1 Introduction 

The physics of cosmic bubble collisions — collisions between bubbles formed from first-order 
phase transitions in models that undergo false-vacuum eternal inflation — has recently at- 
tracted significant attention [1-15]. String theory predicts the existence of a large number 
of meta-stable de Sitter vacua [16-18]. Most of these vacua inflate extremely rapidly, at 
a rate set by the string scale. In a few, the vacuum energy is expected to be comparable 
to the measured dark energy. Hence, in this model the cosmological constant problem is 
solved in the sense that the vacuum energy is an environmental variable that varies from 
place to place in the universe, and in some places it is consistent with the measured value. 

Our observable universe has a very small vacuum energy. In a model that undergoes 
false-vacuum eternal inflation, the most natural origin for such a region is from a bubble 
that nucleated in a first-order phase transition from a (presumably rapidly inflating) parent 
false vacuum. If so, the observable universe fits inside a bubble surrounded by a parent false 
vacuum — a parent vacuum that is evidently unstable to the nucleation of bubbles. Due 
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to this instability, other bubbles are guaranteed to eventually form in the false vacuum 
surrounding ours. Some of these bubbles will collide with ours, making a potentially 
observable imprint on cosmology. 

Much of the recent activity in this area has focused on quantifying the effects of 
collisions on cosmological observables. Significant progress has been made on that front 
[7, 13, 19]. The results obtained thus far took as a starting point a particular form for 
the perturbation at the end of inflation that was first derived in [7] using various approx- 
imations, and the WMAP data set was recently analyzed using the predictions of [7] as a 
template [20]. 

The main result of this paper is a model independent prediction for the collision per- 
turbation in the early universe (Sec. 3). Using the symmetries of the bubble collision as 
a tool, we demonstrate analytically that for a generic perturbation at the beginning of 
inflation that is consistent with the symmetries, the main effect by the end of inflation is 
indeed the one considered in previous work — a curvature perturbation to the early universe 
at reheating that is piecewise linear, but that has a sharp change in its first derivative (a 
"kink") in the direction towards the collision bubble at the edge of the region affected by 
the collision. The kink is sharp in the thin wall limit, and we quantify how it smooths out 
if the collision involves bubbles with thick walls. 

The enormous expansion during the period of slow-roll inflation in the observer's bubble 
after the collision substantially affects the form of the perturbation and gives rise to this 
prediction. By the end of inflation, a "step" discontinuity present at the start of inflation 
evolves into a kink (a discontinuity in the derivative) . It is therefore the sign and magnitude 
of the initial step discontinuity that is the single most important factor in determining the 
cosmological signatures. An analytical approach to this is essential, because — at least for 
the moment — the numerical simulations of bubble collisions that have been done [11, 12, 15] 
cannot incorporate the full period of slow roll inflation. 

The paper is organized as follows. In Sec. 2 we very briefly review the physics of 
bubble collisions, and discuss the approximations and assumptions used in analyzing their 
effects. In Sec. 3 we develop a general technique for analyzing the effects of a bubble 
collision before and during inflation, in Sec. 4 we apply our techniques to several specific 
examples. In Sec. 5 we analyze the effects of collisions involving bubbles with walls of 
finite thickness, and in Sec. 6 we conclude. 

2 Background 

Before taking collisions into account, bubbles that form during a first-order phase transition 
from a false vacuum state will be highly symmetric. The reason for this is twofold: the 
near de Sitter invariance of the false vacuum state, which (since the false vacuum inflates) 
should be very close to perfect, and the fact that the dominant instanton in such a de Sitter 
background is expected to be invariant under rotations in the four Euclidean signature 
dimensions, that is under an 50(4) group of isometries. As a result of this symmetry, each 
unperturbed bubble contains a spatially infinite, homogeneous and isotropic Freidmann- 
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Robertson- Walker (FRW) cosmology with negative spatial curvature (with an SO(3, 1) 
isometry group). Note that none of these conclusions rely on the thin- wall approximation. 

For this cosmology to be consistent with observational data, a period of slow-roll 
inflation is required to solve the flatness problem and produce a scale-invariant spectrum 
of perturbations [1]. The "big bang" of the bubble — FRW time t = 0, where the scale 
factor a(0) = — is a null surface that coincides with the future lightcone of the center 
of the bubble at the moment it appears. If the physics of the phase inside the bubble 
is appropriate, slow-roll inflation takes place along negatively curved time-slices at times 
t > 0; therefore, cosmic bubble collisions occur before slow-roll inflation (Fig. 1). 

A collision between two such bubbles generically makes a large perturbation to the "big 
bang" surface t = 0. Hence, the problem of determining the effects of a single collision on 
cosmology is the problem of evolving an inhomogeneous and anisotropic big bang forward in 
time, which is generally very difficult. Fortunately, the high degree of symmetry preserved 
by the collision is of help. By a de Sitter boost, one can always choose a frame in which 
the two bubbles nucleate simultaneously. In that frame, one expects rotation invariance 
about the axis that connects the centers of the two bubbles, as well as boost invariance in 
the two directions transverse to that axis. In other words, one expects the spacetime and 
field configuration of a 2-bubble collision to preserve an 5*0(2, 1) group of isometries (see 
e.g. [3, 5, 21, 22]). 

In addition, in any model consistent with current observations, the effects of the col- 
lision must be weak and can at some point in the evolution be treated perturbatively. 
Slow-roll inflation must have taken place everywhere in the past lightcone of the earth 
today. As we will see, this significantly simplifies the problem. 

As we will explain, one can choose coordinates in which the collision between the 
two bubbles takes place at one instant everywhere along a hyperboloid. In the thin wall 
approximation discussed below, the effects of the collision are confined to the causal future 
of this hyperbolic surface at that instant. We refer to the (null) boundary of this region as 
the "collision lightsheet." 

2.1 Approximations 

In order to analyze the effects of bubble collisions on observational cosmology, past analyses 
used several approximations and assumptions that we detail here. 

Effective field theory: The collision can be analyzed using a relatively simple effective 
field theory of scalar fields in 3 + 1 spacetime dimensions. This approximation can be 
violated for example in models where the tunneling is from lower dimensions into 3+1 [23— 
25], and its validity in realistic string compactifications has yet to be firmly established. 

Linearity: As we will see, the effects of a single collision (i.e. a collision between two 
bubbles) are fairly easy to determine. Multiple collisions can be analyzed as well provided 
that the interaction between them is weak; that is, provided that by the time the collision 
lightsheets overlap the effects are in the perturbative regime. 
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Dilute bubbles: Prior to any collisions, each bubble is invariant under an SO(3, 1) 
isometry group. The 50(3, 1) invariance holds for bubbles described by the Coleman- 
de Luccia instanton even away from the thin wall limit, as well as for Brown- Teitelboim 
membrane formation [26] and the Hawking-Moss instanton [27]. 

Symmetry breaking: The collision breaks the SO(3, 1) symmetry to an S0(2, 1) sub- 
group (rotation and "boosts" in the plane transverse to the axis connecting the centers 
of the two bubbles). This is the case in simple models for collisions between Coleman-de 
Luccia bubbles, and is supported by numerical simulations [11]. 

Coupling to the inflaton: The collision perturbs the inflaton field, which is a scalar 
or effective scalar mode that drives slow-roll inflation in the bubble after it forms. While 
our analysis does not depend on the specifics of the coupling in any way, it focuses on 
single-field inflation and ignores the effects of the collision on other fields that might be of 
cosmological interest. 

Thin wall: The collision bubble is "thin wall" , meaning that the effects of the collision are 
confined within a sharply defined spacetime region. We will investigate the consequences 
of relaxing this assumption in Section 5. 

Small curvature: In the most convenient set of coordinates, each bubble contains a 
negatively curved FRW cosmology. A single collision produces a "cosmic wake" [19] that 
sweeps across the volume of the bubble. The wake at any given time is bounded by a 
hyperboloid, which previous analyses approximated as planar [7, 19]. As we demonstrate 
in Appendix B, this approximation is valid for our universe, or in general when there has 
been enough slow-roll inflation to solve the flatness problem. 

3 The collision perturbation 

A single bubble, before taking collisions into account, contains a homogeneous and isotropic 
negatively curved FRW cosmology. Before inflation, that is to say at FRW times < 
r < 1/Hi (where Hi is the Hubble parameter during slow- roll inflation) the spacetime 
inside a bubble is dominated by the negative spatial curvature, and can therefore be well- 
approximated as flat Minkowski space sliced with hyperbolic 3-space (eq. (B.2) in the limit 
rHi < 1). For times r > 1/Hi the effects of the inflaton potential become dominant, and 
for the duration of slow-roll inflation the spacetime is de Sitter up to standard slow-roll 
corrections. 

Immediately after a collision, the effects on both the spacetime geometry and the infla- 
ton field may be large. However, if slow-roll inflation begins anywhere in the region affected 
by the collision, the effects of the collision rapidly inflate away and become perturbative. 
Given the lack of large fluctuations in the Cosmic Microwave Background (CMB), we know 
that all the space inside our past lightcone underwent slow-roll inflation. Therefore the 
effects of any collision that is potentially visible today must have become small enough at 
some time to allow inflation to begin everywhere within our past lightcone, and from that 
point in time onward they can be treated perturbatively. This is the key to our analysis. 
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It is important to note that if the inflationary epoch lasted significantly longer than 
the observational bound on curvature requires (number of e-folds N > 62, with the number 
depending logarithmically on various factors such as the reheating temperature) , a bubble 
collision would inflate away into a negligible perturbation in the CMB and would not be 
detected. On the other hand, at least in some simple models [1], the number of e-folds of 
inflation is not expected to greatly exceed the lower bound set by observation. 

We will refer to the time when the effects of the collision become perturbative every- 
where in our past lightcone as to- Knowledge of the inflaton perturbation 5<p and its first 
time derivative at to in this volume provides initial conditions that determine the later time 
evolution of all cosmological observables (see Fig. 1). 

According to our assumption regarding symmetry breaking, the collision perturbation 
preserves an SO{2, 1) isometry group. During slow-roll inflation, the spacetime is approx- 
imately de Sitter. A convenient slicing of de Sitter spacetime that preserves a manifest 
50(2, 1) isometry is: 

-dt 2 

ds 2 = -^y + h(t)dx 2 + t 2 dH%, (3.1) 

where 

h(t) = 1 + (Hit) 2 , (3.2) 

where Hi is the inflationary Hubble scale and dH 2 = dp 2 + sinh 2 pdip 2 is the metric on a 
hyperboloid of unit curvature. 

3.1 Taylor expansion near the lightsheet 

Without a model for the microphysics of the collision, we cannot determine the initial 
perturbation 5(j),5(j) at t = to- Instead, we will simply expand the initial perturbation in 
a general power series and then solve for its time evolution. As we will see, the late-time 
effects are dominated by a single term in the expansion of the initial conditions — the term 
representing a discontinuity in 8(f) across the lightsheet. 

To lowest order in the slow-roll expansion, perturbations 5(p to the inflaton field satisfy 
a free, massless wave equation in de Sitter space [28]: 

U<j> = ~dt [t 2 (l + (m) 2 )^} + 1+ 1 ^0 = 0. (3.3) 

It turns out that this equation can be solved in full generality assuming H2 symmetry, and 
the solution is ([5]): 

6<f>(t,x) = /(£_) - !/'(£_) + <,(£+) + (3.4) 

where / and g are arbitrary functions of one variable, primes denote derivatives with respect 
to the argument, £± = x ± rj are null coordinates, and 

r, = j dt/h(t) = Hr 1 tan _1 (i?^) - vr/(2^) (3.5) 
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Figure 1. Spacetime diagram of a cosmic bubble collision. Each point in the plane of the diagram 
is a two dimensional hyperbolic surface, with a radius of curvature that varies with position in the 
plane. The left-most dashed line bounds the region in the observer's bubble that is affected by the 
collision; referred to in the text as the collision lightsheet. Knowledge of the perturbation 8<j> and 
its derivative provides initial conditions for cosmological evolution in the (red) shaded region. 

is the conformal time. Inflation proceeds for N ~ 60 efolds [29], after which the coordinate 
t ~ H^ 1 e N is very large and rj ~ —(Hft)^ 1 ~ is exponentially small. After a few efolds, 
to a good approximation one can drop the "l"s in (3.2), and (3.4) becomes 

6d>(t, x) « /(£_) + H? V f(£-) + g(£+) - Hfrjg'(U), (3-6) 

where again £± ~ x ± rj, r/ ~ —l/(Hft) is the conformal time. This approximation is 
equivalent to neglecting the effects of spatial curvature on the perturbation and is justified 
in detail in Appendix B. 

To proceed, we wish to describe the perturbation 5(j) in a region near the edge of the 
collision lightsheet, at the time to near the beginning of inflation. As mentioned above, 
in the thin-wall approximation the effects of the collision are entirely confined behind a 
lightsheet: at the initial time t = to, the perturbation 5(p,54> ~ 6(x — xq), where x = xq 



- 6 - 



is the location of the lightsheet at t = to (Fig. 1). (We will explore the consequences of 
relaxing the thin-wall approximation in Section 5.) 

Rather than expanding 5(p, 6(j> at t = to, it is more convenient to expand the functions 
/, g of (3.4). These functions contain the same information as 5<f> and 6<f>. In the thin wall 
approximation they are proportional to a step function that is non-zero only inside the 
collision lightsheet: 

oo 

f(x - r/ ) = M a n H?{x - x ) n 9(x - x ) (3.7) 

n=l 

oo 

g(x + no) =MJ2 b n H?(x - x ) n 6(x - x ). (3.8) 
n=l 

Here a n and b n are dimensionless coefficients, M is a parameter with dimensions of mass, 
and rjo = rj(to) is the conformal time corresponding to t = to- The sum begins with n = 1, 
as an n = term with non-zero coefficient would imply an infinite value for the perturbation 
5(j) (3.6). Written in this form, it is trivial to extend to the full time-dependence: 

oo 

/(x-77) = MY,a n H?(x - x - r] + noTeix - x - v + Vo) (3.9) 

n=l 

00 

g{x + rj) = MJ2 b n H?(x -x + rj- Vo) n 0(x - x + r? - r/ ). (3.10) 

n=l 

These expansions determine 5(j)(rj,x). After many efolds of inflation when rj ~ 0, 
g is non-zero everywhere inside the lightsheet, x — xo > f]o = —\vo\- By contrast, the 
right-moving mode / is non-zero only in the region of spacetime x — xo > —r]o = |?7o|> a 
coordinate distance Ax = 2 1 770 1 horn the edge of the lightsheet. Physically, / corresponds 
to a component of the perturbation 5(f) that moves directly away from the observer at the 
speed of light, starting from early in inflation (illustrated in Fig. 1 by the dashed line in 
the £;+ direction originating from the point (xo,to))- As such it is unobservable, at least in 
a situation in which the collision lightsheet intersects the observer's last scattering volume. 
Therefore, in computing the observable effects of a bubble collision, we can neglect the 
rightmoving modes /(£-). 

As can be seen from (3.9), in the vicinity of the lightsheet at late times (x ~ xo+rjo, rj ~ 
0) the dominant term in (3.9) is the one proportional to b\ (barring a large hierarchy 
among the dimensionless coefficients b n ). Longer inflation suppresses the higher terms 
exponentially, since it reduces the range Ax within the observer's last scattering volume 
by a factor that scales with e~ N (it also suppresses the effects of the term proportional to 
b\ — but by the lowest power if e~ N ). 

Therefore, generically the late-time perturbation is 

5c/)(ri,x) = MHibix6{x) + 0{rj,x 2 ), (3.11) 

where for notational convenience we have chosen the origin of x so that xo = —rjo and the 
collision lightsheet is at x = when rj = 0. 
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To relate this to the perturbation at early times rj = rjo, note that = b\Hi£ i+ 9(£ + ) 

corresponds to a perturbation 6(j)(x,rjo) = MHib\x6(x + rjo). This contributes a disconti- 
nuity 



in the inflaton field value at the edge of the collision lightsheet. Hence, the amplitude 
of the "left-moving part" of the discontinuity in 8(f) at the edge of the lightsheet is what 
determines the coefficient of the x6{x) kink at late times (this same conclusion is derived 
in a slightly more direct way in Appendix A) . The parameter Mb\ , along with the angular 
size of the affected disk and its position on the CMB sky, fully parametrize the leading 
effects of a thin-wall bubble collision. 

Our conclusion is that a general initial inflaton perturbation evolves to a perturbation 
that has been smoothed essentially by integrating in x once (see Appendix A). In other 
words, an initial discontinuity becomes a kink, a kink becomes a function with discontinuous 
second derivative, etc. This is in accord intuitively with the more standard treatment of 
inflationary perturbations in fe-space, where 5(f) ~ k~ 3 / 2 (l ±irjk). Near the end of inflation 
\r]k\ <C 1, and so an initial perturbation in Fourier space evolves to one with one less power 
of k. 

4 Examples 

As we have seen, the dominant effect of a bubble collision arises when the collision creates 
a discontinuity in the inflaton field across the collision lightsheet (when 61 in (3.8) is not 
zero). When present near the beginning of slow-roll inflation, such a discontinuity in 
the field "inflates away" into a discontinuity in the first derivative of the field (x6(x)), 
which was the starting point for various analyses of the effects on cosmological observables 



In this section we discuss two approaches in which we can analytically study the 
effects of the collision on the universe before inflation. While either approach could in 
certain microphysical models be used to determine the coefficients a n and b n in the initial 
perturbation to the inflaton, doing so is not our primary purpose. Instead, we wish to show 
that the effect of a bubble collision on the inflaton across a broad class of models is in fact to 
create a discontinuity in its value across the lighsheet (in the thin wall approximation) — in 
other words, we wish to show that b\ 7^ generically, rather than calculate its value in any 
specific model. 

Our first approach makes use of the "free passage" approximation of [12], while the 
second treats the domain wall between the bubbles as a boundary condition on the inflaton. 
In a sense, these two approaches constitute opposite extremes — free passage applies when 
the kinetic energies are so large that the fields can be treated as free, and the domain 
wall between the bubbles is irrelevant or even non-existent, while treating the domain 
wall as a boundary condition requires strong non-linearities in the field. Nevertheless, 
in both cases, the collision generically produces a discontinuity in the inflaton across the 
collision lightsheet. In addition, numerical simulations support the conclusion that thin- 



A(f> = MHibxlrjo 




[5, 13, 19, 30]. 
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wall bubble collisions lead to discontinuities in the field value along the lightsheet of the 
collision [9, 11, 12, 15]. 

4.1 Free passage 

When two planar scalar-field domain walls collide head-on, if the kinetic energy of the 
walls exceeds the relevant energy differences in the scalar potential, the field configuration 
immediately after the collision is determined in a very simple way by the pre-collision field 
values [12] (see Appendix A of [31] for an extension to collapsing spherical configurations 
that "turn inside out"). Consider a model with several scalar fields and a potential V((f>). 
Initially, a region of = 0o separates two regions where = 0i and = 02 bounded by 
ingoing planar walls (see Fig. 2). The results of [12] show that immediately after the walls 
collide, the field value in region 3 is 

03 = 01 + 02 - 00- (4.1) 

In the case of free fields (V(4>) constant) in flat Minkowski spacetime, the field in region 
3 will remain fixed and equal to 03 for all times after the collision. The same follows in 
interacting models if 03 happens to be a minimum of the potential. In general, the post- 
collision evolution of the field in region 3 depends on the model — given = 03 as an initial 
condition, the time evolution of 0(t) will depend on the scalar potential ^(0), as well as 
the spacetime curvature. 

We wish to consider a collision between two thin-wall bubbles containing field values 
01 and 02, embedded in a rapidly inflating parent vacuum 0o- As before we will ignore the 
hyperbolic curvature of the walls (although the analysis of [31] can be extended to that 
case). Bubble 1 is the observation bubble, and we will assume 0i falls on or in the basin 
of attraction of a slow-roll inflationary slope that ends with an observationally acceptable 
vacuum. Bubble 2 is the collision bubble, and 02 may be a minimum or close to a minimum 
of the potential. 

We would like to apply this approach to a model that is potentially consistent with 
observation. This requires the observer's bubble to undergo a period of slow-roll inflation. 
A simple such model is the single-field model illustrated in Fig. 2, where the observer's 
minimum = O & S is separated from the parent vacuum by a long slow-roll plateau. After 
tunneling, the fields near the bubble walls are not in the minimum. Instead, the tunneling 
point — which is the value of the field on the walls of the bubble — is separated in field space 
from the false vacuum only by the width of the barrier. Cosmological evolution inside 
the bubble is rapid at first, but "friction" due to spatial curvature prevents "overshoot" 
and causes the inflaton to move only a distance roughly equal to the barrier width before 
coming to near rest and then slowly rolling along the plateau [1]. 

If the kinetic energy in the walls is large at the time of collision, the free passage result 
(4.1) should apply immediately after the collision. The field value = 03 may or may not 
fall on or (in a multi-field model) in the basin of attraction of the slow-roll inflationary 
plateau of bubble 1. If it does not, the collision will have a very large effect on the region 
3, which is almost certainly not compatible with observation. Therefore we will assume 
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that it does fall on the plateau, so that standard slow roll inflation takes place inside the 
affected region. 

In this case the fields in regions 1 and 3 may behave approximately as free fields for 
at least a time of order the inflationary time determined by the height of the slow- 
roll plateau, at which point slow- roll inflation begins. The analysis of Sec. 3.1 applies 
directly to this situation, with the initial field values approximately determined by (4.1). 
We assume that the backreaction of the energy in field gradients is small, so that the 
spacetime in regions 1 and 3 is de Sitter space with Hubble constant H^. 

The simple equation (4.1) must be applied with care to collisions involving such bub- 
bles, because the field inside the bubble is evolving rather than constant. However, for 
early times after the collision and for the first few efolds of inflation it should be valid in 
the vicinity of the tunneling point. The most appropriate value for the free-passage approx- 
imation is one not far from the tunneling point (labelled c\ and C2 in Fig. 2). Numerical 
simulations are probably required to determine the accuracy of this approximation in any 
given model. 

With 0o = O,0i = ci, and 02 = C2 (see Fig. 2), the value of the field in the region 
affected by the collision (region 3 in Fig. 2), is (4.1): 

03 - ci + c 2 (4.2) 

soon after the collision. The evolution of the field in region 2 is unobservable, but for 
illustrative purposes in this section we will treat region 2 as de Sitter space with the same 
vacuum energy as regions 1 and 3. 

To find the time-dependence in region 3, we will make use of the left and right moving 
formalism developed in the first section. Expanding as in (3.4), we should require that 
/ and g, are non-zero only in the regions 1 and 2 respectively, and that both are non-zero 
in region 3. Additionally, they should satisfy 

/(*) " h'M = C2 (4.3) 
g(z) + ^g'(z) = ci (4.4) 

in the regions where are non-zero; (x c ,t c ) are the coordinates of the collision point. 
The solutions to these equations are 

f(z) = c 2 (l-e t ° z ) (4.5) 
g(z) = Cl (l-e-^ z ) (4.6) 

and therefore the solution for with the complete causal structure is 

4>(x, t) = c 2 (l- e^-^+'Jc) + k e tc(x-v-xc+r,c)\ (_ x + V + Xc _ v j + 
ci ^1 - e-Mz+^c-^) + h e -t c (x+ v -x c - Vc )^ e ^_ x _ v + Xc + rj( j 

(4.7) 
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Figure 2. Left panel: A potential V{<j>) with three minima. The Coleman-de Luccia tunneling 
paths from the central false vacuum are indicated with arrows. Bubbles that tunnel to the right 
will inflate along the slow-roll plateau leading down into the observation bubble's minimum. Right 
panel: Schematic of the field configuration prior to the collision. 

where 77 = — 1/t and r/ c = — l/i c . It is straightforward to check that this solution shows 
the correct behavior around the collision point (x c ,t c )- 

At late times t > 1 (corresponding to many efolds of inflation), the last term in the 
parenthesis of each line of (4.7) becomes negligible, leaving only the (1 — e"'). The leading 
term of this expression is linear in x, showing again that the late time behavior near the 
light cone has a discontinuity in the first derivative, but not in the field value. 

In Fig. 3 we show an example of a collision between two identical bubbles in the 
free-passage approximation. The figures in table 3 gives snapshots at different times of the 
field in (4.7). 

4.2 Predictions for the temperature of the spot 

At least in models involving scalar fields, even minimal knowledge of the field space and the 
free-passage approximation can lead to a prediction for the sign (hot or cold) of the spot in 
the CMB temperature map. In the case of a collision between two identical bubbles, (4.1) 
shows that the initial perturbation always pushes the field away from the tunneling point. 
For example, in a collision between two bubbles of the right-hand minimum in Fig. 2, the 
analysis above suggests that the field value in region 3 will be fa ~ 2c\. 

In a model where fa ~ 2ci > (p b s (4>obs is the minimum in the observer's bubble) the 
consequences for cosmology in region 3 are likely to be large and are unpredictable without 
more detailed knowledge of the model — the field could swing back across the minimum 
4> bs and onto the inflationary plateau, or it could immediately reheat in the minimum. In 
any case it seems unlikely that the affected region (region 3) will undergo more efolds of 
inflation than the unperturbed region (region 1). However, if the plateau is long compared 
to the barrier width, the field in region 3 may advance relative to the field in region 1 by 
the equivalent of only a few efolds — enough to have a potentially observable effect, but 
not necessarily a very large one. The result will be fewer efolds of inflation in the affected 
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Figure 3. The figures are snapshots of a symmetric collision at different times. In this set up, 
Ci = C2 = 1, while the false vacuum is at = 0. The collision time is at to = 15, and the figures show 
snapshots at times just before the collision (ti = 8), just after (f 2 = 20), a while after (i 3 = 100) 
and a long time after the collision (£4 = 1000). We notice that the perturbation freezes (due to de 
Sitter background) and that the jump in <f> disappears at late time leaving only a kink, a jump in 
the first derivative. Notice that the four figures are not at the same scale. 

region, since the perturbation pushes the inflaton forward down the inflationary plateau. 
In both cases the result is fewer efolds of inflation in the affected region. After these results 
were obtained and while the manuscript was in preparation, [15] and [32] appeared that 
perform numerical analyses in specific models that appear to support this conclusion. 

Fewer efolds of inflation produces an intrinsic cold spot, heuristically because such a 
region reheats earlier and has longer to cool before last scattering. Due to the Sachs- Wolfe 
effect this produces a hot spot on the CMB (and a corresponding azimuthal polarization 
pattern around the edge of the disk [13]). This conclusion is not model-independent — for in- 
stance, non-linear effects from interactions could invalidate the free-passage approximation — 
but we expect it to hold in a fairly wide range of models. 

By contrast, in a collision with a bubble containing another type of vacuum, the jump 
in the inflation can have the opposite sign. A realistic model of the string landscape would 
contain a multi-dimensional field space. Assuming ^3 is in the basin of attraction of the 
slow-roll inflationary trajectory, the effects are (tentatively) more likely to delay inflation 
in the affected region for the same reasons explained above. In principle, they could also 
lead to interesting features in the power spectrum of perturbations in the affected region 
of the CMB sky — such as non-Gaussianities or oscillations — which would be produced as 
the field oscillates around the "valley floor" inflationary trajectory (a somewhat related 
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phenomenon was recently studied in e.g. [33]). 
4.3 Domain wall boundary conditions 

An alternative approach is the one followed in [7], where the domain wall separating the 
two bubbles after the collision was treated as a boundary condition for the inflaton. The 
simplest boundary conditions are Dirichlet — that the inflaton field takes a fixed value on the 
domain wall. Given the trajectory of the domain wall (which can be determined by treating 
it as a thin wall separating two regions of de Sitter space [5, 34]), this fully determines 
the evolution of the inflaton field everywhere inside the observation bubble. The initial 
conditions for the inflaton prior to the collision are set by the Coleman-de Luccia initial 
conditions along the lightcone of the nucleation point of the observer's bubble, perturbed 
by the presence of the domain wall (Fig. 1). Together, the Coleman-de Luccia initial 
conditions plus the domain wall form a complete Cauchy surface for the evolution of the 
inflaton. 

The inflaton perturbation during slow-roll inflation is determined by the functions 
/(£_) and <?(£+) of (3.6). Consider solving the boundary value problem described in the 
paragraph above. The inflaton perturbation 5(f) due to the collision is zero outside the 
collision lightsheet. Because the coordinate £_ is constant along null rays crossing the 
collision lightsheet, this means the function / must be zero or constant. Therefore, the 
inflaton perturbation is determined by the single function <?(£+). 

Another way to understand this argument is to note that the value of 5(f) at any point 
in the £+, plane is determined by its values at any two points, one along each of the past 
directed null geodesies originating from the point. Since the left-going past null geodesic 
terminates on the Coleman-de Luccia initial condition surface, it is only the right-directed 
past null geodesic that determines 6(f), and this in turn determines the function <7(£+). 

The above arguments demonstrate that in this scenario, the perturbation depends only 
on the coordinate £ + perpendicular to the lightsheet of the collision. Since the perturbation 
is zero outside the lightsheet, in the thin wall approximation the function <?(£+) must be 
proportional to — £o)> as i n (3.6). The question of whether or not 5(f) is actually 
discontinuous — rather than simply having discontinuous first or higher derivatives — comes 
down to precisely what boundary condition one imposes at the domain wall. While we will 
not present further details here, in various examples it is clear that for generic choices, 5(f) is 
itself discontinuous. Put another way, because the initial conditions form a Cauchy surface, 
requiring 5<p to be continuous on the lightsheet overconstrains the problem. Therefore, 
generically there is a non-zero discontinuity in 6(f) (which then inflates into a discontinuity 
in the first derivative as explained in detail in Sec. 3 and App. A). 

Our analysis on this point differs from that of [7], which required continuity of the 
perturbation along the lightsheet. However, our conclusion — that the inflaton perturba- 
tion due to a thin-wall collision near the end of inflation has a discontinuity in the first 
derivative — is identical. The sign of the effect in the case of Dirichlet boundary conditions 
depends on the value of the inflaton on the domain wall in a way that was discussed in [7] . 
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5 Thickness of the wall 



In this section we will explore the consequences of relaxing the thin-wall approximation. 
Specifically, we will study a scenario where the collision bubble has walls of finite thickness. 
In such a case the perturbation 5<j) 18 no longer proportional to theta function, but instead 
to a smooth function related to the field profile of the collision bubble's wall when it formed 
in the parent false vacuum. 

In our analysis we treat the thick wall as if it were composed of two thin walls separated 
by a distance w when the collision bubble forms in the parent false vacuum, and determine 
the length scale d corresponding to w in the inflaton perturbation at the beginning of 
inflation (Fig. 4). 

We consider two distinct limits. Choosing a reference frame in which the two bubbles 
nucleate simultaneously and the observer's bubble is centered at x = 0, we first study the 
case where the collision bubble nucleates close to ours (at a distance Hjp < 1); then we 
consider the scenario in which the colliding bubble nucleates close to the de Sitter horizon 
of the false vacuum (n — Hjp <C 1). 

In order for the collision to be visible, the lightsheet of the collision should intersect 
the last scattering surface of the observer. That means that the lightsheet must intersect 
the FRW timeslice at the beginning of inflation near the comoving location of the observer. 
As can be seen from Figures 1 and 4, when the collision occurs early as it does in our first 
limit (p small), the lightsheet intersects the the beginning of inflation surface far to the 
left of x = 0, and therefore an observer that can see it is boosted away from the collision 
bubble. By contrast, if the collision bubble nucleates sufficiently close to the horizon of 
the parent false vacuum, the collision lightsheet may intersect the beginning of inflation 
surface near x = and hence could be visible to an observer at rest in this frame. 

5.1 Notation and setup 

We work with conventions illustrated in Fig. 4. Due to the symmetry of the Coleman-de 
Luccia initial conditions, prior to the collision all parts of the wall — defined for instance 
by contours of constant field value — will follow a hyperbolic (constant proper acceleration) 
trajectory in the parent false vacuum (curved red line segments in Fig. 4). Approximating 
the observer's bubble by the lightcone of its nucleation point for the moment, a and b mark 
the events where the inner and outer surfaces of the collision bubble's wall collide with 
ours. 

After the collision, we propagate the effects along future directed null lightsheets (the 
straight red segments in Fig. 4). These lightsheets intersect the beginning of inflation 
surface at the points marked a and f3, separated by a distance d. Our goal in this section 
is to determine d from w, the collision bubble's radius r, and its center location p. Because 
the spacetime curvature of our bubble is negligible until the beginning of inflation at a time 
of order H~ l 3> HJ 1 , we treat the interior of our bubble as Minkowski spacetime. 
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Figure 4. Colliding with a thick-wall bubble centered at x — p. The red lines indicate the inner 
and outer edges of the thick wall and the region it affects. As usual, each point in this diagram 
represents a 2-hyperboloid with a radius that depends on position in the plane of the diagram. 



5.2 Collision soon after nucleation 

In this first limit, the two bubbles are nucleated at a distance p -C HJ (Hubble distance 
of the false vacuum), the initial size of the bubble is r c and the thickness of the wall at 
nucleation is w. In this limit, as explained above we can ignore the curvatures of both 
the parent false vacuum and slow roll inflation and treat the spacetime globally as flat 
Minkowski. 

We set the coordinates in such a way that our bubble nucleates at x = 0. Then the 
right branch of light cone of our bubble is given by x = t, and therefore the equations to 
solve to find x a and x b are 

-x 2 a + (x a -p) 2 = (r c + w) 2 (5.1) 
-x 2 + (x b - p) 2 = r 2 c , (5.2) 

where we have accounted for the fact the the wall follows a hyperbolic trajectory. Even if 
at leading order we have x a ~ Xb — HJ 1 /2, the difference is 

x a ~ Xb — • (5.3) 

p 
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We compute the coordinates of the point a (f3 would be analogous) by finding the 
intersection between the hyperbolic surface at the end of inflation and the causal cone of 
a; namely we have to solve the system 

- x 2 a + t 2 a = H7 2 (5.4) 
x a 2x a = t a . (5-5) 



We find 



H~ 2 - Ax 2 

x a = (5-6) 

Ff- 2 _|_ A T 2 



and the result for (xp,tp) would be the same with a o b. Therefore 

X a -X/3 = (x a - X b ) ^1 



H7 2 



H7 2 



(5.8) 



t a -tp = (x a -x b )[l + ^- Xb ). (5.9) 



The (hyperbolic) distance is now easily found: 

d 2 = (Ax) 2 - (At) 2 



= H7 2 {Xa Xb) \ (5.10) 



X a Xl) 



which can be written 



<?~ H - 2 (2r c w + w 2 ) 2 

- i ;p 4_ ;p 2( r 2 + ( rc + u; )2) +r 2( rc + u; )2- ) 

In order to compare with the computation we will show in the next section, it is useful 
to introduce the dimensionless parameter 5 = Hjp (the distance between the bubble centers 
in units of HJ 1 ): 

d2 = S (2r cW + W 2 ) 2 

Hf 5 A -5 2 H 2 f (rl + (r c + w) 2 ) + Hjrl(r c + w) 2 ' 1 ' } 

For our analysis to be valid the bubbles should not overlap when they nucleate, meaning 
that p > r c + w. Taking the limit p S> r c ,w gives 

dHi = H 2 f w(2r c + w)5' 2 (6 < 1). (5.13) 

The quantity dHi is the apparent wall thickness in units of the inflationary Hubble length; 
as such it roughly corresponds to the angular size of the apparent thickness on the CMB 
sky after minimal inflation. 
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5.3 Collision with a bubble nucleated close to the de Sitter horizon 

When p <C HJ 1 the effects of the parent vacuum's curvature prior to the collision are 
negligible. By contrast, when the bubble nucleates close to the horizon of the de Sitter 
false vacuum (ttHJ 1 ) the effects of the de Sitter curvature become important. 

Again making use of the hyperbolic metric (3.1), the trajectories of the wall and the 
light cone of our bubble respectively are given by: 

cos z (x — (7T — ej j cos z {x + e) 
x = tan" 1 ^), (5.15) 

where for simplicity we set Hf = 1 (we will restore Hf restore at the end) and e = ^HJ 1 —p 
is the distance from the de Sitter horizon where the bubble nucleates. Also, from now on 
we might use c(r c ) to indicate the constant cos 2 (r c ). 

Equating (5.14) and (5.15) to find the collision point gives 

1 c(r c ) 

o7 \" = oT \ * (0.10) 

cos z [x) cos A (x + e) 

The point x = tt/2 is midway between the center of the observer's bubble and the de Sitter 
horizon, and the divergence there indicates that a collision at that point would take place 
only after infinite time. Since we are interested in collisions with bubbles that nucleate 
near the horizon, we expand in7 = 7r/2 — x 1 and e <C 1. One finds 

This leads to only one physical solution once we impose the condition that 7 — > when 
r c — > 0, so we find 

1 — cos(r c ) 
sin 2 (r c ) 



7 = ^E^ (5-18) 



and therefore 



7r 1 — cos(r c ) 
2" sin 2 (r c ) 

1 — cos(r c ) 
sin 2 (r c ) 



x& = — —y- — (5.19) 



h = ( - ZT V > ) (5-20) 



and analogously one can find 



7r 1 — cos(r c + w) 



(5.21) 



t a = coi(- \.,;y c Z: J e). (5.22) 



2 sin 2 (r c + u>) 

1 — cos(r c + w) 
sin 2 (r c + w) 

We need now to match these results with the coordinates inside our bubble, that we 
treat as flat, and then determine the intersection points of the lightsheets with the surface 
at the beginning of inflation. Inside our bubble we use the hyperbolic slicing of Minkowski 

ds' 2 = -dt' 2 + dx' 2 + t' 2 dH' 2 2 (5.23) 
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so that the matching condition with the outside metric (3.1) is that the hyperboloids need to 
have the same radius on the light cone, where the two patches meet: t'\i c = t\i c . Therefore 

< = t' a = t a (5.24) 

x b = t' b = tf,. 

It is easy now to follow the null geodesies from these points and find their intersection 
with the beginning of inflation surface: 

_ t' 2 + x' 2 = Hr 2 (5.25) 
x' 2 - 2t a = -t' (5.26) 



giving 



TT-2 

*a = ^+ta (5-27) 

4 = "^+ta (5.28) 



and 



jj-2 

X ' p = ^- + t b (5.29) 

rr-2 

fy = —±- + t h . (5.30) 

Computing the hyperbolic distance on the surface at the beginning of inflation is 
straightforward: 

2 (ta - tb) 2 



d = H ~i 

tah 

r4 



H} (2r c w + w 2 ) 2 (1 + C ot(H f e/2)) 2 2 2 
H 2 64 cot 2 (if /e /2) /£ 



Hf (2r c w + w 2 ) 2 2 2 
H 2 64 f 



H 2 f e 2 (5.31) 



where we expanded for small r c and w (r c ,w ^ H ^ 1 ) and in the last line we expanded for 
small Hj 1 e as well. 



Introducing again the parameter 5 = Hjp the equation (5.31) can be rewritten 

dHi~ H 2 f w(2r c + w) 7 ^- (tt-<5<C1). (5.32) 



Comparing with (5.13) reveals the effects of the observer's boost. When the bubbles 
nucleate very close together, the observer must be moving away from the collision bubble 
with a large Lorentz factor in order to be in position for the collision lightsheet to intersect 
her last scattering volume. This boost away from the collision has the effect of Lorentz 
dilating the apparent thickness of the wall, which accounts for the negative power of 5. By 
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contrast, if the collision bubble nucleates very close to the de Sitter horizon, the observer 
must be boosted towards it, and the thinner it appears. 

Most bubble collisions will occur with a 8 > 1. For such a case the apparent thickness 

is 

dHi ~ H 2 f w(2r c + w)/8 (6 ~ 1). (5.33) 

For a thin wall bubble Htw -C 1, and in all cases Hf{r c + w) < 1 (the critical size of the 
bubble must be less than one Hubble length). 

The results of [19] showed that a thin-wall bubble produces a very interesting feature 
on the CMB sky at the degree scale (a "double peak" in polarization intensity, for in- 
stance). The quantity dHi is the ratio of the apparent width of the wall to the radius of 
curvature of the universe at the initial time. Ratios of length scales are fixed during FRW 
expansion; therefore the right-hand side of (5.33) is also the ratio of the wall thickness to 
the curvature at the time of last scattering. Observational constraints on curvature show 
that < .01 in the current epoch [35]. Using fij. = (aH)~ 2 and values from WMAP-7 
[35], the observational constraint combined with equation (5.33) shows that 

d L s/r sh > 300H 2 f w(2r c + w)/8, (5.34) 

where d^s is the lengthscale corresponding to the wall thickness at last scattering, and 
r s h is the sound horizon. The double peaks of [19] are separated by an angular distance 
of order twice the sound horizon. Therefore in order for the thin-wall approximation to 
accurately describe such a feature, the left-hand side of (5.34) should be less than one, and 
therefore 

{H f w)H f {2r c + w)< 1(T 2 . (5.35) 

For collision bubbles with thicker walls, the thin wall approximation does not suffice to 
describe degree-scale features. For instance, the double peak is likely to be strongly affected 
or smoothed over entirely. Given the profile of the collision bubble's wall (which could be 
computed from the underlying microphysical model), a detailed analysis could be done to 
determine precisely what the effects on temperature and polarization are near the edge of 
the affected disk. If there is substructure within the wall, it could have an interesting and 
non-trivial effect on cosmological observables in the region "inside" the wall. Conversely, 
the detection of a bubble collision in the CMB would allow one to infer something about 
the characteristics of the walls of the collision bubble, a remarkable fact considering its 
truly microscopic origin. 

6 Conclusion 

The detection of a cosmic bubble collision would be a discovery of tremendous importance. 
As such, it is important to quantify its signatures accurately. Because of the lack of a 
unique or easily calculable microphysical model, it is difficult to make entirely unambiguous 
predictions. Nevertheless, we have shown here that the primary effects of bubble collisions 
in a broad class of models are parametrized by only a few parameters, and that this 
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conclusion is robust under relaxing at least some of the assumptions made in previous 
analyses. 

One of the remarkable aspects of this problem is the potential for discovering certain 
details of near-Planck length physics by observing enormous structures that extend across 
the entire observable universe. One example is the color (hot or cold) of the affected disk 
in the CMB temperature map, which we have argued is indicative of the type of vacuum in 
the collision bubble. Another example is the thickness of the collision bubble's wall, which 
could imprint in an interesting way on the cosmological signature. 

One major direction that remains largely unexplored is to study a bubble collision 
in string theory, or in a model with at least some of its ingredients (compactified extra 
dimensions, fluxes, charged branes, etc.). One obstacle is that a microphysical model 
for slow-roll inflation would be required, since the interactions of the inflaton field are of 
major importance in determining the signatures. We believe that the analysis of Section 
3 is general enough to cover most such cases, at least when slow-roll inflation is effectively 
single-field. However, in multi-field or more exotic inflation models interesting features 
might emerge. Another potential signature is the effect of the collision on other light 
fields, which could potentially lead to perturbations in physical "constants" such as the 
fine-structure constant a. or probes of early universe physics such as 21cm [36]. Finally, 
recent work on "barnacles" is interesting and should be further explored [37]. 
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A Series expansions 

In the analysis of Sec. 3, we expanded the functions / and g in a power series in x at an 
early time rj = r/o during inflation (Eqs. 3.7 and 3.8). The coefficients in this expansion 
determine the initial perturbation 5<j) at rj = rjo, and the full ij dependence of /, g and 5<p 
is trivial to determine from the coefficients on and bi given (3.4). However, under some 
circumstances it is more natural to instead expand the perturbation 5(p at r/ = r/o in a 
power series (for instance, if the initial perturbation is determined by the microphysics). 
In this appendix, we determine the evolution of 5(j) in rj in terms of the coefficients in its 
expansion at rjo. 

To begin, by analogy with a single term in (3.8), consider a perturbation that is a 
monomial in x and zero outside the collision lightsheet: 



Because only the left-moving part of the peturbation is observable, for simplicity we will 
set /(£-) = (this condition implicitly determines the initial time derivative d v 5(p). Under 
these conditions, to find the evolution we wish to solve 
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The general solution, in terms of an integration constant, reads 

C 



9 P (x + Vo) = e^VWrft (Y(l + p) - T (l + p; 



Vo 



e (*-*°)/«>T$(r{i+p)-r(i+p, 



X — Xo 

Vo 



9(x - x Q ). 



(A.3) 



Here x) = t s ~ 1 e~ t dt is the incomplete gamma function. 

There are two conditions (A.3) must satisfy. When we take the derivative, the term 
with the 5 function arising from the differentiation of the Heaviside 9 function must vanish: 
this is the case for all values of C. Second, the term which does not multiply the Heaviside 
9 function should vanish identically, because this solution should be non-zero only in the 
region affected by the collision (for x + rj > xq + t]o); this can be accomplished by choosing 
the integration constant C = 0. 

With C = 0, the solution for all r\ is 

J_ (i+ - xo - Vo 
ml 



g p (H + ) = (-f +1 p\ e «+-«o)/.o _ £ 



m=0 



e^-xo-w), (A A) 



where as before = x + V- For notational simplicity, we can (by a shift in the coordinate 
x) choose xq = — 770, and will do so for the remainder of this appendix. 

Consider the lowest order term p = 0, which is the solution for the initial perturbation 
<j)(x — xq) = cq9{x — xq) for some constant cq. The solution for g is simply 

9o (£+) = co(l - < 

and the field perturbation is 

<%))(£+) = co(! " 
For late times {rj — > 0) this becomes 

5<p m (^)^c (l-e^)9(U)- (A.7) 

As we found in section 3, an initial discontinuity in 5(f) inflates into a continuous 
function with discontinuous first derivative at the collision lightsheet. 
In general, the initial condition at rjo can be expanded 



Vo 



(A.5) 
(A.6) 



5(f>(x, vo) = ^2 Cn ( x ~ x o) n ®( x ~ x o) 

n=0 

By linearity, the solution at general rj is then 

oo 



(Ai 



Safari) = Y,cn(-) n+1 n\ 



n=0 



J+/V0 



E 1 



m=0 



ml \ r/o 



V_ 
Vo 



E 

m=0 



[m — 1)1 



+1 n! 



n=0 



n 

'°-E i 

Z J Til 



m=0 



(m-iy 



ml \ r/o 



0(6 



(A.9) 



where the last line is valid as ry — > 0. 
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B Spatial curvature 

As mentioned in Sec. 2.1, the curvature of the collision lightsheet becomes small in after 
inflation. In fact, we will see that its radius of curvature is always larger than the radius 
of curvature of the FRW slice. 

Using the coordinates (3.1) the radius of curvature of the lightsheet is simply t. The 
radius of curvature of the FRW slice at the beginning of inflation is most easily found using 
a different set of coordinates: 



cosher) = y/l + (Hit) 2 cos(Hix), (B.l) 
sinh(iJjr) sinh /3 cos 9 = a/TT {Hit) 2 sin(Hix), 

smh(HiT) cosh (3 = Hit cosh p, 
sinh(iTjT) sinh (3 sin 9 = H{t sinh p, 
<P = f- 

In these coordinates the metric (3.1) becomes 

ds 2 = -dr 2 + Hr 2 sinh 2 HirdH$, (B.2) 
dflf = d(3 2 + sinh 2 f3(d9 2 + sin 2 Odip 2 ). (B.3) 

In these coordinates, after N efolds inflation ends homogeneously at time r e ~ NH^ 1 , and 
the radius of curvature of the universe at that time is sinh(i7jT e ). We wish to compare this 
to the curvature of the lightsheet at the same time. To do this, we should first find the 
intersection of the collision lightsheet — which is well approximated by the lightcone of a 
bubble that nucleated at x = xq — with the hypersurface r = r e . Using (B.l): 

V 7 ! + (Hit) 2 cos(^x) = cosh^Te) ~ € -^- (B.4) 
y/1 + (Hit) 2 cos(Hi(x -x )) = l (B.5) 

where we used the approximation T e Hi = N S> 1. 

Taking the ratio of the two equations above one gets 



cos( Hix) e HiTe 



COs(Hi(x - Xq)) 



(B.6) 



Since t > r S> H i 1 for all x (c.f. (B.l)), equation (B.5) gives t ~ H i 1 / cos(Hi(x — xq)) 
and therefore 

2 cos(Hix) 

showing that the radius of curvature of the collision lightsheet at the end of inflation is of 
order or larger than the radius of the universe at that time. The denominator cos(Hix) 
goes to zero when x ~ ±H~ 1 ir/2, which happens when the collision bubble nucleated very 
close to the horizon of the false vacuum (+) or very close to the observer's bubble at x = 
(")■ 
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